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The notion of "paired" fermions is central to important condensed matter phenomena such as su- 
perconductivity and superfluidity. While the concept is widely used and its physical meaning is clear 
there exists no systematic and mathematical theory of pairing which would allow to unambiguously 
characterize and systematically detect paired states. We propose a definition of pairing and develop 
methods for its detection and quantification applicable to current experimental setups. Pairing is 
shown to be a quantum correlation different from entanglement, giving further understanding in the 
structure of highly correlated quantum systems. In addition, we will show the resource character of 
paired states for precision metrology, proving that the BCS states allow phase measurements at the 
Heisenberg limit. 



I. INTRODUCTION 

The notion of pairing in fermionic systems is at least 
as old as the seminal work of Bardeen, Cooper and Schri- 
effer explaining superconductivity The formation of 
fermionic pairs with opposite spin and momentum is not 
only the source for the vanishing resistance in solid state 
systems, but it can also explain many other interesting 
phenomena, like superfluidity in helium-3 or inside a neu- 
tron star. 

For instance, with recent progress in the field of ultra- 
cold quantum gases fermionicpairing has gained again a 
lot of attention @, i, i, 1, 1, H S ■ These experiments 
allow an excellent control over many parameters inher- 
ent to the system, offering a unique testing ground for 
existing theories and an exploration of new and exotic 
phases. However, the notion of pairing in these systems 
is less clear and sometimes even controversial. Recent 
experiments on the BEC-BCS crossover have caused a 
heated debate whether or not the obtained data was in 
agreement with pairing I?! [ToL [ill . [l2| . In addition, pair- 
ing without superfluidity [131 ] has been observed in these 
experiments, raising fundamental questions on quantum 
correlations in fermionic many-body systems. 

Motivated by these exciting experiments and the cen- 
tral role pairing plays in many physical phenomena, and 
by the perceived lack of accepted criteria to verify the 
presence of pairing in a quantum state, we propose a 
clear and unambiguous definition of pairing intended to 
capture its two-particle nature and to allow a system- 
atic study of the set of paired states and its properties. 
We employ methods and tools from quantum informa- 
tion theory to gain a better understanding of the set of 
fermionic states that display pairing. In particular, we 
develop tools for the systematic detection and for the 
quantification of pairing, which are applicable to cur- 
rent experiments. Our approach is inspired by concepts 
and methods from entanglement theory, thus building a 
bridge between quantum information science and con- 
densed matter physics. 

Since they contain non-trivial quantum correlations, 



paired states belong to the set of entangled many-body 
states. However, pairing will turn out to be not equiva- 
lent to any known concept of entanglement in systems of 
indistinguishable particlesfli, El, EJ, El EE El EE HI, 
HI EE El EE EE EE EE Ejbut to represent a particular 
type of quantum correlation of its own. We will show that 
these correlations can be exploited for quantum phase es- 
timation. Hence pairing constitutes a resource in state 
estimation using fermions as much as entangled states 
with spins. 



This article is organized as follows: After the intro- 
duction of the language necessary for the description of 
fermionic systems in Sec. UH we will introduce the gen- 
eral framework of pairing theory in Sec. IIIII This part 
includes our definition of pairing and methods for its de- 
tection and quantification. In order to fill the theory with 
life we will apply it to two different classes of fermionic 
states in Sees. IIVI and [Vj We start out with pairing in 
fermionic Gaussian states in Sec. IIVI The interest in 
this family of states is two-fold. First, the pairing prob- 
lem can be solved completely in this case, so that Gaus- 
sian states are particularly interesting from a conceptual 
point of view. Second, there exists a relation between 
pure fermionic Gaussian states and the BCS states of su- 
perconductivity (see Sec. Ill Dl for the details) which are 
examples of paired states par excellence. This enables 
us to translate methods developed for the detection and 
quantification of pairing for Gaussian states to the BCS- 
states. The reader interested in the application of our 
pairing theory to experimental application is referred to 
Sec. El There we study pairing for number conserving 
states, i.e. states commuting with the number operator. 
This class includes the states appearing in the BEC-BCS 
crossover, and we will develop tools for the detection of 
pairing tailored for these systems. In I VII we will show 
that certain classes of paired states constitute a resource 
for quantum phase estimation, proving that pairing is a 
resource similar to entanglement. 
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II. FERMIONIC STATES 

In this chapter we review the basic concepts needed 
for the understanding of fermionic systems. We start out 
with some notation used for the description of fermionic 
systems in second quantization in Sec. Ill Al As pairing 
is a special sort of correlation, we continue with a review 
on quantum correlations and entanglement in systems 
of indistinguishable particles in Sec. Ill Bl This general 
part is followed by the introduction of fermionic Gaussian 
states and number conserving states in Sees. Ill CI and 
III Dl The latter includes the introduction of BCS-states 
and their relation to the Gaussian states. As this part is 
only necessary for the application of the pairing theory to 
these concrete examples in Sees. IIVI and IVl it is possible 
to skip this part at the beginning, and then refer to it 
later on. 



A. Basic notation 

We consider fermions on an A/-dimensional single 
particle Hubert space TL — (D M . All observables are 
generated by the creation and annihilation operators 
a] and aj,j = 1,...,M, which satisfy the canonical 
anti-commutation relations (CAR) {a,k,ai} = and 
{a,k,a\} — 5ki- We say a] creates a particle in mode 
(or single particle state) ej, where {e^} C <D M denotes 
the canonical orthonormal basis of TL. In general, for any 
normalized / £ TL, we define af = ^ fc fj(ij, the annihi- 
lation operator for mode /. 

Sometimes a description using the 2M hcrmitian Ma- 



form in the complex representation a& 



J2i Ukiai, 



jorana operators c^j-i = a] + a,j, C2j = (— *)( a j — a j)i 
which satisfy {cfc,q} = 25m, is more convenient. 

The Hilbert space of the many-body system, the an- 
tisymmetric Fock space over M modes, Am, is spanned 
by the orthonormal Fock basis defined by 



\ni,...,n M ) 



(«!)' 



|0>, 



(1) 



where the vacuum state |0) fulfills aj |0) = Vj. The rij £ 
{0, 1} are the eigenvalues of the the mode occupation 
number operators rij = a\a,j. The A^-particle subspace 
spanned by vectors of the form ([TJ) satisfying ^\ rij = N 
is denoted by . The set of density operators on the 

Hilbert space TL = Am,A^ is denoted by S(TL). 

Linear transformations of the fermionic operators 
which preserve the CAR are called canonical transforma- 
tions. They are of the form Ck >— > c' k = J^OfcjCj, where 
O £ 0(2M) is an element of the real orthogonal group. 
These transformations can be implemented by unitary 
operations Uo on Am which are (for detO = 1) gener- 
ated by quadratic Hamiltonians in the Cj (see, e.g. [30|). 
The subclass of canonical operations which commute 
with the total particle number A^p = m are called 
passive transformations. They take a particularly simple 



where U is unitary on the single-particle Hilbert space TL, 
i.e., they describe (quasi)free time-evolution of indepen- 
dent particles. Canonical transformations which do not 
commute with N op are called active. They mix creation 
and annihilation operators. 



B. Quantum correlations of fermionic states 

The notion of " pairing" used in the description of su- 
perconducting solids, superfiuid liquids, baryons in nu- 
clei, etc. is always associated with a correlated fermionic 
system. The sub ject of quantum correlations in fermionic 
systems is vast [31(. In recent years, there has been re- 
newed interest from the perspective of quantum informa- 
tion theory. There quantum correlations (aka entangle- 
ment) of distinguishable systems (qubits) play a crucial 
role as a resource enabling certain state transformations 
or information processing tasks. The detailed quantita- 
tive analysis of quantum correlations motivated by this 
has proven to be valuable also in the understanding of 
condensed matter systems (see [32| for a review). 

In contrast to the usual quantum information setting 
which studies the entanglement of distinguishable par- 
ticles, the indistinguishable nature of the fermions is of 
utmost importance in the settings of our interest. The ex- 
isting concepts for categorizing entanglement in systems 
of indistinguishable particles fall into two big classes: En- 
tanglement of modes [13, EE EE E3, EE Ea, 123, H3, E2I 

and entanglement of particles. Entanglement of particles 



2EJ2EJ23, m, m 



has been considered e.g. in [23|, |24 
leading to the concept of Slater rank [24|, [25| , being the 
generalization of the Schmidt rank to indistinguishable 
particles. We show in Sec. IHII that our definition of pair- 
ing does not coincide with any of the existing ideas. We 
refrain from giving an exhaustive review on the existing 
concepts referring the interested reader to the mentioned 
literature and references therein, and restrict to the fol- 
lowing definition: 



Definition II. 1. A pure fermionic state p p 



(N) 



l^p V ^)( x ^p V ''l e £(-4m ') * s called a product state, if there 
exists a passive transformation a k 



a' k such that 



N 



i«#°> = IKi°>- 

i=i 



(2) 



A state p s is called separable, if it can be written as the 
convex combination of product states, i.e. 



K 



o - V A o (AW 



(3) 



0=1 



£ S(A { m p) ) are 



where J2 P =i A p = 1, A p > and all p p c ^v-^Af 
product states. All other states are said to have "Slater 
number larger than 1" and are called entangled (in the 
sense of Jm, \Mj). 
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We denote the set of all separable states by S sep and 



by S. 



(N) 



of particle number N . 



S sep H S(Am ) the set of all separable states 



Note that the sets S sep ,Ssep of separable states are 
convex and invariant under passive transformations. 
Both properties will be useful later on. 
Separable states have only correlations resulting from 
their anti-symmetric nature and classical correlations due 
to mixing. In the terminology of Refs. 0, [H| they have 
Slater number one and describe unentangled particles. 
These states will certainly not contain correlations as- 
sociated with pairing. (Note that they can be mode- 
entangled for an appropriate partition of modes.) 

Besides basis change, there are other operations, which 
do not create quantum correlations and it is useful to see 
that the set of separable states is invariant under them. 

Lemma II. 2. Let p G S sep be a separable state. Then 
the state after measuring the particle number nh — a \ a h 
in some mode h is separable for both possible outcomes 
n h =0,1. 

Furthermore, ph = tr ah [p] , the reduced state obtained by 
tracing out the mode ah, is also separable. 

Proof. As S sep is convex, it is sufficient to prove the claim 
for product states p. Let \*S?) = Yi^Li a / ■ |0) be the vector 
in Hilbert space corresponding to p. Our aim is to show 
that |\E r ) = |\E'o) + |5'i), where \^i) are product states and 
n/i = I eigenstates of the occupation number operator 
nh- If h is in the span of {/i<fc<jv} or orthogonal to 
it, the state already is a nh eigenstate and we are done. 
Otherwise, define /jv+i orthogonal to the fk<N such that 
h G span{/i<fc<j\r +1 } and define another orthonormal 
basis {gj} for the span with g\ = h and 52 oc /j\r+i — 
(h ■ fN+i)h [here (h ■ /jv+i) denotes the inner product 
on the single particle Hilbert space]. Then we can write 

i*> = °/* + i°!w + i n £i4,i > = ^ +^) n f= + iHi°> 

for some x,y € C. Hence, |*) = |* ) + with |*o) = 
xUf+^allO) and = -ya^V^ |0) which both 

clearly are product states. The reduced state tr/ t [|^)(^ r |] 
is the statistical mixture of |^o) and and therefore 
clearly separable. □ 



C. Fermionic Gaussian states 

Fermionic Gaussian states are represented by density 
operators that are exponentials of a quadratic form in 
the Majorana operators. A general multi-mode Gaussian 
state is of the form 



p = K exp 



c T Gc 



(4) 



anti-symmetric matrix can be brought to a block diagonal 
form 



M 



OGO- = ( ° "A 



(5) 



by a special orthogonal matrix O G SO(2M). 

From eq. ([4]) it is clear that Gaussian states have an in- 
terpretation as thermal (Gibbs) states corresponding to 
a Hamiltonian H that is a quadratic form in the Cfe, i.e., 
H = jC T Gc = jJ2k>i G ki[ck,ci] and the form eq. © 
shows that every Gaussian state has a normal-mode de- 
composition in terms of M single-mode "thermal states" 
of the form ~ exp(— (3a*a). From this one can see that 
the state is fully determined by the expectation values of 
quadratic operators a^a.,- and a\aj. These are collected 
in a convenient form in the real and anti-symmetric co- 
variance matrix Y which is defined via 



T ki = -tr (p[c k ,ci\). 



(6) 



It can be brought into block diagonal form by a canonical 
transformation: 



M 



OYO T = 



Xj 
-A, 



(7) 



For every valid density operator, Xj G [—1,1], and the 
eigenvalues of Y are given by HXj. Hence, every Y cor- 
responding to a physical state has to fulfill iY < 1 or, 
equivalently, rr^ < 1 and to each such Y corresponds 
a valid Gaussian density operator where the relation be- 
tween G and Y is given by Xj = tanh(/?j/2). The co- 
variance matrix of the ground state of H is obtained in 



the limit \@j 



00 i.e., A,- 



sign(/3 3 ). In fact, this 



shows that every pure Gaussian state is the ground state 
to some quadratic Hamiltonian. The purity of the state 
can be easily determined from the covariance matrix as 
a Gaussian state is pure if and only if T 2 = — 1 (see. e.g., 
0). 

As mentioned, Gaussian states are fully characterized by 
their covariance matrix and all higher correlations can be 
obtained from Y by Wick's theorem (see, e.g., [HI]) via 



i p tr[pcj 1 



L '32p 



] = PffowJ 



(8) 
is the 

^2 _ 



where 1 < j\ < . . . < j% p < 2M and r jl: ... j2p 
corresponding 2p x 2p submatrix of Y. Pf(r j ' 1) ... i3 - 2 ■ — 
det(r jli ... j2p ) is called the Pfaffian. 

In some cases it is more appropriate to use a different 
ordering of the Majorana operators, the so-called q-p- 
ordering c = (c x , c 3 , . . . , c 2 m-V,c 2 , c 4 , • . • , c 2 m), opposed 
to the mode-ordering introduced at the beginning. When 
using the q-p-ordering, the relation between the real and 
complex representation is given by 



where c = (ci, . . . , Qjjvf), if is a normalization constant 
and G is a real anti-symmetric 2M x 2M matrix. Every 



2 T = n a T , 



n = 



1 

it 



1 

-it 



(9) 
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,t 



,t 



). The transformation 



where a = {a\ &m, a[ , . 
matrix ft fulfills ftftt = 21. 
In the q-p-ordering the covariance matrix obtains the fol- 
lowing block structure: 



r 



(10) 



Finally, for some purposes it is more convenient to use 
the complex representation, where the covariance matrix 
is of the form 



1 



r — -O^fn — . 

,. - LU y R Q 



Q R 



(11) 



where Q k i = (i/2[a k , a/]), R k i = {i/2[a k ,a}]) and Q de- 
notes the complex conjugate. Note that R^ = —R and 
Q T = -Q and hence = -T c . The condition ff t < 1 
takes the form 4r c rJ < 1. 

The description of p by its covariance matrix is espe- 
cially convenient to describe the effect of canonical trans- 
formations, i.e. time evolutions generated by quadratic 
Hamiltonians: if c k i— > OkiQ in the Heisenberg picture 
then r i— > OTO T in the Schrodinger picture. For a pas- 



sive transformation a k 



J2i Ukl a li the q-p-ordered 



Majorana operators transform as 



JT 



O n c 1 



0„ 



X Y 

-Y X 



(12) 



where X = Re(£7) is the real part of the unitary U, and 
Y = lm(U) the imaginary part. Note that O p is both 
orthogonal and symplectic. The behaviour of T c under 
a passive transformation is particularly simple: Q and R 
transform according toQ« UQU T and R t— > URW . 

Passive transformations can be used to transform pure 
fermionic states to a simple standard form, the so-called 
Bloch-Messiah reduction [H]. The q-p ordered CM f B cs 
takes the form (JTUJ) where 



r _ _ r _ (T\ ( -2Im(u fc w*) 



(13) 



r - f ~ \v k \ 2 2Re(u k v* k ) \ 
1 9P " W{-2Re(u k v* k ) \u k \ 2 -\v k f ) 



In Hilbert space, the state in standard form is given by 



I* 



(N) >. 
Gauss' 



H(u k +v k atal k )\0), 



(15) 



where u k ,v k £ €, \u k \ 2 + \v k \ 2 = 1, N = Y, k ( a l a k) = 
2 J2 k \ v k\ 2 - This comprises the kind of "paired" states ap- 
pearing in the BCS theory of superconductivity [l| with 
k = (k, t), — k = (— k, I). We will refer to these states 
as Gaussian BCS states. We would like to stress the fact 
that every pure Gaussian state is a Gaussian BCS state 
in some basis. 



D. Number conserving fermionic states 

For the application to physical systems we are inter- 
ested in states for which the particle number is a con- 
served quantity. We call p a number conserving state if 
[p, N op ] = where -/V op denotes the total number opera- 
tor. Thus, the density operator of a number conserving 
state can be written as a mixture of A^ op -eigenstates. In 
particular, all separable states as defined in Def. Ill, ll are 
number conserving. 

The Gaussian BCS wave function (| 1 5[) is not number 
conserving (except for the case J2 k \ u kV k \ — that either 
u k or v k vanishes for every mode) , but a relation to these 
states can be established via the identity 



2M 



(16) 



N=0 



where the number conserving 2iV-particle BCS state is 
given by 



\<%)=C N ^ ak P^ 



N 



(17) 



where we have introduced the pair creation operator 
Pfc = a\a}_ k . The coefficients a k are related to u k and v k 
via oifc = v k /u k , Cn is a normalization constant which is 
seen to be 



-1/2 



Cn = (AH) 2 J2 



\a jl | 2 . . . \cti |i 



j !<...< j N 



by rewriting Eq. (|17[) as 



C N N\ 



jl<h<— <3N 



a kl 



*k N Pl...PL\0). (18) 



The coefficients Xn — ([\ k u k ) /(N\Cn) can be inter- 
preted as the probability amplitude of being in state 

I^bcs) smce J2n I Aiv 1 2 = 1 - We will in general drop the 
term number conserving and refer to states of the form 
(fl?) as BCS states. 

Whenever the distribution of the An is sharply peaked 
around some average particle number TV, expectation val- 
ues of relevant observables for the number conserving 



BCS states I* 



(AT) 

BCS, 



are approximated well by the expec- 



tation values of the Gaussian BCS state. This relation 
will turn out very useful later on, as results on Gaussian 
states can be translated into results on number conserv- 
ing BCS states. 



III. PAIRING THEORY 

In this section we introduce a precise definition of pair- 
ing as a property of quantum states. 
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A. Motivation and statement of the definition 

The simplest system in which we can find pairing con- 
sists of two particles and four modes 1 . The prototypical 
paired state, for example the spin-singlet of two electrons 
with opposing momenta, is of the form 

l$) = i(44+44)|o). (19) 

The states describing many Cooper pairs in BCS theory 
are generalizations of |$). 

The state |$) describes correlations between the two 
particles that cannot be reproduced by any uncorrelated 
state and it can be completely characterized by one- and 
two-particle expectations (consisting of no more than two 
creation and annihilation operators each) . This is a char- 
acteristic of the two-particle property "pairing" that we 
propose to make the central defining property of paired 
states in the general case of many modes, many parti- 
cles and mixed states. Since, moreover, we would call 
the state |$) paired no matter what basis the mode op- 
erators a, refer to and we want it to comprise all BCS 
states, we are led to the following list of requirements 
that a sensible definition of pairing should fulfill: 

1. States that have no internal quantum correlation 
must be unpaired. These are the separable states 
©• 

2. Pairing must reveal itself by properties related to 
one-and two-particle expectations only. 

3. Pairing should be a basis-independent property. 

4. The standard "paired" states appearing in the de- 
scription of solid state and condensed matter sys- 
tems, i.e., the BCS states with wave function (TT71) 
must be captured by our definition. 

Further, it would be desirable that there exist examples 
of paired states that are a resource for some quantum 
information application. 

Let us define: 

Definition III.l. The set of all operators {O a } a on Am 
which are the product of at most two creation and two 
annihilation operators is called the set of two-particle 
operators. We denote it by A^. 

These operators capture all one- and two-particle prop- 
erties of a state p and should therefore contain all infor- 
mation about pairing. We will call a state p paired, if it 
can be distinguished from separable states by looking at 
observables in A^ alone. This is formalized in the follow- 
ing definition: 



For three modes, all pure two-particle states are of product form. 



Definition III. 2. A fermionic state p is called paired 
if there exists a set of operators {O a } a C A2 such that 
the expectation values {tr[pO a ]} cannot be reproduced by 
any separable state p s € S sep . 
States that are not paired are called unpaired. 

This definition automatically fulfills our first two re- 
quirements by definition. The third, basis independence, 
clearly holds, since the set of separable states is invariant 
under passive transformations. We will show that the 
last requirement is met, both for Gaussian and number 
conserving BCS states, i.e. all of them are paired (see 
Lemma IV. 31 and Subsec. HVBp . Moreover, in Sec lVIl 
we can show that there exist paired states that are a 
resource for quantum metrology. 

For states with a fixed particle number, i.e. p G 
S(Afy[ ) it is sufficient to compare with expectation val- 
ues on iV-particle separable states pa € Ssep , as for 
all other states the expectation values of rij) and 
((Si n «) 2 ) differ due to the particle number constraint. 
For number conserving states, only number conserving 
observables lead to non- vanishing expectation values and 
one can thus restrict to linear combinations of 4 a i' 

44 a * a J- 

For Gaussian states pairing must reveal itself by prop- 
erties of the covariance matrix, as all higher correlations 
can be obtained from it via Eq. ([5]). This important 
fact enables us to give a complete solution of the pairing 
problem for fermionic Gaussian states, which we present 
in Sec. |IVl 



B. Relation of pairing and entanglement 

Paired states are fermionic states exhibiting non-trivial 
quantum correlations. In particular, by definition paired 
states are inseparable i.e., entangled in the sense of 0, 
[25j . This raises immediately the question: Is pairing 
equivalent to entanglement? Below, we provide examples 
of entangled, but unpaired states that demonstrate that 
pairing is not equivalent to entanglement (of particles) 
but represents a special type of quantum correlation 2 . 

Lemma III. 3. There exist states that are entangled ac- 
cording to the Slater rank concept, but not paired. 

Proof. Consider the state | M/4 ) = |(444 a l + 
4444)|0), which is entangled according to the Slater 
rank definition. However, one sees immediately that 



2 Note that our basis-independent definition clearly has no relation 
to entanglement of modes, which is basis-dependent. The prod- 
uct states of Def. Ill.ll can be mode-entangled for some choice of 
partition of modes, e.g. l/v / 2(l| + a|)|0) is entangled in modes 
4 and a\. 
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the one-and two-particle expectations for |\&4) are the 

same as for pi 4) = + 5 1 ^> 2 > (^2 1 , where |$i) = 

al^Oga^O), ($2) = a i a e a 7 a s\^)- Since p^ is a product 
state, ISI/4) is not paired. One can construct further ex- 
amples in a similar manner using e.g. other states with 
higher Slater rank. □ 

Since pairing is defined via expectation values of 
one-and two-particle operators only, one might wonder 
whether pairing is related to entanglement of the two- 
particle reduced state. To study this relation, we recall 
the definition of the two-particle density operator and the 
closely related two-particle density matrix (see eg. [35j): 

Definition III. 4. Let p be the density operator of a 
fermionic state. Then O^-w^n = tr[pa\a^aiak] is called 
the two-particle reduced density matrix (RDM). It is 
usually not normalized and fulfills tr[0^} — (N^ ) — 
(N op ). 

The operator P2 — /tr[0^] is called reduced two- 
particle density operator (RDO). 

Note the crucial difference between the two-particle 
RDM and the RDO. While the RDM contains all two- 
particle correlations of p, the RDO corresponds to the 
two-particle state of any two particles when the rest of 
the system is discarded. We would like to emphasize, that 
pairing is not equivalent to entanglement of the RDO, 
and therefore it is a property of the one- and two-particle 
expectations: 

Lemma III. 5. Let \^gcs) ^ e a num ber conserving BCS 
state as defined in (|17p with ak = 1 V k = 1, . . . , M . 
Then its two-particle RDO (see Def. \III-4\ ) PescS2 * s 

always paired. However, p^g(j S 2 * s entangled if and only 
if M :i .V 2. 

The proof is given in Appendix lEl 

We would like to stress the point that Lemma IIII.5I 
shows the existence of paired states that are not entan- 
gled. Having assured that our definition of pairing does 
not coincide with entanglement, we now turn to methods 
of detecting and quantifying pairing. 



C. Methods for detecting pairing 

Taking Def. IIII.21 we aim at finding tools that can 
be used for the detection and quantification of pairing. 
These will be applied to systems of Gaussian states and 
number conserving states in Sees. IIVI and fVl respectively. 
In this section, we exploit the convexity of the set of un- 
paired states to introduce witness operators and obtain a 
geometrical picture of the set. The quantification of pair- 
ing via pairing measures will be discussed in Sec. IIIID1 

Given a fermionic density operator, we are interested 
in an operational method to determine whether it is 



paired or not. As in the case of separability, this simple- 
sounding question will turn out to be rather difficult to 
answer in general. 

Starting from Def. IIII.21 it is clear that the set of 
unpaired states is convex. This suggests the use of the 
Hahn-Banach separation theorem as a means to certify 
that a given density operator is not in the set of paired 
states. In analogy to the entanglement witnesses in quan- 
tum information theory [3(| we define 

Definition III. 6. A pairing witness W is a Hermitian 
operator that fulfills tr[Wp u ] > for all unpaired states 
p u , and for which there exists a paired state p such that 
tr[W p] < 0. We then say that W detects the paired state 
P- 

The witness defines a hyperplane in the space of den- 
sity operators such that the convex set of unpaired states 
lies wholly on that side of the plane characterized by 
triply] > 0. According to the Hahn-Banach theorem 
[37| , for every unpaired state there exists a witness oper- 
ator which detects it. In principle, a witness operator can 
be an operator involving an arbitrary number of creation 
and annihilation operators. However, since definition of 
pairing refers only to expectation values of operators in 
A2, it is enough to restrict to operators from that set. 
This represents a significant simplification both math- 
ematically (witness operators from a finite dimensional 
set) and experimentally, since operators involving more 
than two-body correlations are typically very difficult to 
measure. 

The construction of entanglement witnesses detecting all 
entangled states is an unsolved problem in entanglement 
theory, and we will not be able to give a complete solu- 
tion to the problem of finding all pairing witnesses either. 
However, in Section [V] we will construct witnesses for a 
large subclass of BCS-states by using the correspondence 
betweenF number conserving and Gaussian BCS states. 

Whether a state p is paired or not can be determined 
from a finite set of real numbers, namely the expectation 
values of a hermitian basis {O a } of A^. This allows us to 
reformulate the pairing problem as a geometric question 
on convex sets in finite-dimensional Euclidean space, de- 
scribe a complete set of pairing witnesses, and deduce a 
relation to the ground state energies of quadratic Hamil- 
tonians. 

Consider a set {O a , a = 1, . . . , K} C A2 of hermitian 
operators in A2 that are not necessarily a basis. Denote 
by O the vector with components O a . We define the set 
of all expectation values of O for separable states 

C d = {v = tv{6p s ] : Ps £ S sep } C R A '. (20) 

For a state p let v p = tr[Op]. By definition, p is paired 
if v p £ Cq. As the set of separable states is convex, so 
is Cq . Hence, we can use a result of convex analysis to 
check if v p ^ Cq (see e.g. [Hj]): 
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Lemma III. 7. Let C C R be a closed convex set, and 
letve R N . Then 

v E C O Vr e WL N : v- r > E(r) = inf w ■ f. (21) 

For our purposes, this translates in the following result: 

Lemma III. 8. For a vector of observables O — 
(0 1 ...,0 K ) let H{r) = r ■ 6 and E{r) = 
mi peSscp {tr[H(r)p}}. Then W(f) = H(f) - E(r) is a 
pairing witness, whenever E(f) ^ iu.f a \i p {tr[pH(f)]}. 

If{O a } form a basis of A2, then W(r) is a complete set 
of witnesses in the sense that all paired states are detected 
by some W{f), i.e., p is unpaired iff tr[W(r)p] > OVf. 

Proof. The witness property of W(r) is obvious from the 
definition of E(r). 

For the second part, "if" is clear and "only if" is seen 
as follows: By Lemma |III.7[ if tr [W(r)/o] > OVr then 
v p G C, i.e. the expectation values can be reproduced 
by a separable state. But since all expectation values of 
operators £ A2 can be computed from v p this implies all 
two-particle expectations of p can be thus reproduced, 
i.e. p is unpaired. □ 

For an Af-mode system with annihilation op- 
erators a,-, a standard choice of O a is, e.g., 
given by the real and imaginary parts of 
{{a\a]a k ai) i>0 . k> i,{a\a]) l>j ,{a\a j )i>j}, i.e., the 
dimension of A2 (as a real vector space) is 
K = M 2 (M - l ) 2 /2 + 2M 2 . 

Thus Lemma IIII.8I gives a necessary and sufficient cri- 
terion of pairing and provides a geometrical picture of the 
pairing problem. While the proof that a state is unpaired 
will in general be difficult as it requires knowledge of all 
E(r) and experimentally the measurement of a complete 
set of observables, practical sufficient conditions for pair- 
ing can be obtained by restricting to a subset O C Ay,. 
We will show in Sec. IV Al that for a certain choice of 
{O a } C Ay the set Cq has a very simple form and allows 
a good visualization of the geometry of paired states and 
the detection of all BCS states up to passive transforma- 
tions. 

To provide a way to determite E(f) used in 
Lemma IIII.8I we point out an interesting connection to 
the covariance matrices T c (cf. Eq. (fTTj) ) of Gaussian 
states: even for number conserving states, E(f) is given 
by a quadratic minimization problem in terms of r c . 

Lemma III. 9. Let E(r) and H(r) be as in Lemma [LLL.8\ 

and let O — {ajajafcd;, a\aj} and group the components 
of r in two subsets (r*)ijki and (r)ij corresponding to 
the one- and two-particle observables, respectively. Then 
E{f) is given by a quadratic minimization problem over 
complex covariance matrices Eq. ill]) , in particular the 
off-diagonal block R ofT c . We have 

E{r) = inf {7 T Af(r)7 + w{rf^} , (22) 

R= — J?t 

4fl 2 = -l 



where (7)^ = (aj.a ; ) = -iRik + j5ki and the r- 
dependent quantities are [M (r)](ik)(ji) = ^^ijki + r^ijik 
and [w(f)]ki = r kl . 

The minimization can be extended over all ( not necessar- 
ily pure separable) CMs without changing the result. 

Proof. The minimum rmn P £s si!p {(H (r)) p } is attained for 
pure separable states, i.e., product states. All pure 
fermionic product states are Gaussian; then by Wick's 
theorem the expectation values of the Oijki — a\a -akO,i 
factorize as (ala^a k ai) p = (ajat)(ajtaj) - (a\a k ) (ata fe ) + 

(a|ct;) (Ojdk) ■ Since product states are also number con- 
serving, the first term vanishes. For the other two we use 
that (atai) = —iRik + j$ku i-e., they only depend on the 
off-diagonal block R. The pure state condition T 2 = — 1 
translates into 4i? 2 = — 1 for product states Q = 0. 
We could extend over all CMs 7 C since only the block R 
appears in the expression to be minimized over and since 
if r c (Q, R) is a valid CM then so is T c (0, R). □ 

This lemma provides a systematic way to construct 
pairing witnesses. 



D. Pairing measures 

It would be desirable if a theory of pairing not only 
answers the question whether a state is paired or not, 
but also quantifies the amount of pairing inherent in a 
state. For this purpose, we introduce the notion of a 
pairing measure: 

Definition III. 10. Let p be an M -mode fermionic state. 
A pairing measure is a map 

M: p^ M{p) 6 JR+, 

which is invariant under passive transformations and ful- 
fills Ai{p) = for every unpaired state p. 

In addition, it is often useful to normalize A4 such 
that A4(po) = 1 defines the "unit of pairing". The pair 
state of Eq. (fT9|) would be an obvious choice for this 
unit, but as we see in Sec. IIVDI for Gaussian states a 
different unit is more natural, therefore we do not include 
normalization in the above definition. 

In the geometric picture of the previous section, a can- 
didate for a pairing measure that immediately comes to 
mind is the distance of v p from the set C. This measure 
is positive, and it is invariant under passive transforma- 
tions, as those correspond to a basis change in the space 
of expectation vectors. 

The computation of this distance is, in general, very 
difficult and there is no evident operational meaning to 
this quantity. In the following sections we will introduce a 
different measure that can be computed for relevant fam- 
ilies of states and allow a physical interpretation in terms 
of quantifying a resource for precision measurements. 
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IV. PAIRING FOR GAUSSIAN STATES 

In this section we study pairing of fermionic Gaussian 
states. We start with the construction of pairing wit- 
nesses in Sec. IIV Al which will later be a useful guide- 
line for the construction of pairing witnesses for number 
conserving states. Then we derive a simple necessary 
and sufficient criterion for pairing of Gaussian states. In 
Sec. IIV CI we show how pure fermionic Gaussian states 
can be connected to an SU(2) angular momentum repre- 
sentation. This picture will guide us to the construction 
of a pairing measure. 

A. Pairing witnesses for Gaussian states 

Pairing witnesses for pure Gaussian states emerge nat- 
urally from the property that every such state is the 
ground state of a quadratic Hamiltonian (see Sec. Ill Cp . 
This leads to the following theorem: 

Theorem IV. 1. Let < e < 1 and let < \v k \ 2 < 1 - e 
and \ vk\ 2 > 0. Then the operator 

M 

H = ^2(l-e- M 2 )(n fe +n_ fc ) 
fe=l 

- 2v k u* k Pl ~ 2v* k u k P k (23) 
is a pairing witness, detecting 

l^Gauss) =Y[{u k +V k Pl)\0). 

k 

Proof. Every Gaussian state is the ground state of a 
quadratic Hamiltonian. In particular, \^ Gauss) is seen 
to be the ground state of 

M 

Ho = ^2(\u k \ 2 - \v k \ 2 )(n k +n_ fe - 1) 
fc=l 

- 2v k u* k Pl - 2v* k u k P k 

with the help of (fTB|) - (fTi|) . as the Hamiltonian matrix 
of Hq and T can be brought simultaneously to the stan- 
dard forms |[SJ) resp. ([7]). Subtracting the minimal en- 
ergy for separable states, E^f n = -(1 - 2e)J2 k ( n k + 
n- k ) — (\u k \ 2 — \v k \ 2 ) (note that (P k ) = for separable 
states), we arrive at the Hamiltonian ([23]) . For sepa- 
rable states p, the expectation values of P k vanish, so 
that tr[Hp] > 0. For the Gaussian BCS state, however, 
{^Gauss\H\^ G auss) = -4e£ fe \v k \ 2 < 0. □ 

B. Complete solution of the pairing problem for 
fermionic Gaussian states 

Every Gaussian state is completely characterized by 
its covariance matrix, so that the solution of the pairing 
problem must be related to it. The pairing problem is 
completely solved by the following theorem: 



1 



l 


unpaired states 






* J „ J_„ J \ / 


7 ► 







FIG. 1: Bloch sphere representation of the expectation values 

of j k x \ an d 3k^ f° r a variational BCS state. All pure 
states lie on the surface of the sphere. Unpaired states lie 
on the z-axis, while the maximally paired states lie on the 
equator. 



Theorem IV. 2. Let p be the density operator of a 
fermionic Gaussian state with covariance matrix T c de- 
fined in (jlip . Then p is paired iff Q ^0. 

Proof. First, note that the condition Q — is indepen- 
dent of the choice of basis. If p is not paired, then there 
exists a separable state having the same covariance ma- 
trix as p. This implies Q = 0, as separable states are 
convex combinations of states with fixed particle num- 
ber, and thus (i/2[a k , at]) = 0. 

Now, let r c be the covariance matrix of a paired Gaus- 
sian state, and assume that Q = 0. As R is anti- 
hcrmitian, there exists a passive transformation such that 
Rij = TiSij, and Q = is unchanged. But such a covari- 
ance matrix can be realized by a separable state fulfilling 
(ni) — ri in contradiction to the assumption. □ 



Note that Thm. IIV. 2 1 implies that a Gaussian state is 
unpaired iff it is number conserving. 



C. Angular momentum algebra for Gaussian states 

In this section we will show that pairing of Gaussian 
states can be understood in terms of an SU(2) angular 
momentum algebra. The expectation values of the angu- 
lar momentum operators can be visualized using a Bloch 
sphere, giving us further understanding of the structure 
of pairing in Gaussian states. It later leads to the con- 
struction of a pairing measure for these states. Define 
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the operators [39l. |4C 



Jk 



Av) _ 

Jk — 

A') _ 

Jk — 



n- k ) 



They fulfill j { k a ',j { k ' = ie a bcj k , a,b,c S {x,y,z}, 
forming an SU (2) angular momentum algebra. For pure 
Gaussian states in the standard form (fT"5| the expec- 
tation values of the angular momentum operators are 
given by (jf>) = Re(u k v k ), {j { k y) ) = lm(u k v k ), and 

Of) = |(1 - 2\v k \ 2 ). As j 2 = Ei^yjjf) 2 = 1/4 
independent of u k and v k , the expectation values for ev- 
ery pure Gaussian state lie on the surface of a sphere 
with radius 1/2. As we have shown in Thm. IIV.2[ ev- 
ery unpaired state p u fulfills (j k x ^) Pu = (j k )p« = ^' so 
that these states are located on the z-axis. The states 
on the equator have (j£ ) 2 + (j k ) 2 = 1/4, i.e. they 
correspond to \u k \ 2 — \v k \ 2 = 1/2. The situation is de- 
picted in Fig. [TJ Referring to the states on the equa- 
tor as maximally paired is suggested by the fact that 
they have maximal distance from the set of separable 
states. This intuitive picture is further borne out by 
two observations: first, the states on the equator dis- 
play maximal entanglement between the involved modes 
[30]. Second, they have the property 3 that they achieve 
the minimal expectation value of any quadratic witness 
operator (up to basis change). To see this, recall from 
Sec. Ill CI that any quadratic Hamiltonian of two modes 
k, —k is (up to a common factor and basis change) of 
the form at + s'm6(n k + fi- k ) + cos#(Pj + P k ). It is 
a witness (i.e., has positive expectation for all product 
states), if a > |max{0, 2 sin6*}| and does detect some 
paired state as long as sinf? > — 1. The minimum eigen- 
value is sin# — 1 + a and the minimum tr (Wp) = — 1 is 
attained for p = |(1 + P^)|0)(0|(1 + P k ). 

The pairing measure which is the topic of the next 
section will confirm the characterization as maximally 
paired. 



Definition IV. 3. Let p be a fermionic state, and let 
Qki = i/2tr(p[a k , ai]). Then we define 



M G (p) 



kl 



(24) 



Lemma IV. 4. j\4g as defined in Eq. (IV. 3\) is a pair- 
ing measure fulfilling J\4g(p) < M for every M-mode 
Gaussian state. 

Proof. Under a passive transformation Q i— > UQU T , and 
hence 1 1 <5 1 1 2 i s invariant. Further, we know by Thm. IIV.2I 
that Q — for unpaired states. 

It remains to show that for an M-mode Gaussian state p 
we have Mc(p) < M. Let T c be the 2M x 2M covariance 
matrix of p defined in (fTTjl . Wc show first that M.(p) is 
maximized for pure Gaussian states. To do so, recall that 
an admissible covariance matrix for a Gaussian state in 
the real representation has to fulfill iT < 1 with equality 
iff r is the covariance matrix of a pure Gaussian state. 
This translates into r c rj < 1 with equality iff r c belongs 
to a pure Gaussian state. Using the form of F c given in 
(fTTj) , this implies 2tr[QQt + RR^} < tr[l] = 2M. Hence, 
I Mil < M — \\R\\ 2 . It follows that for fixed value of 
\\R\\ 2 the value of ||Q||| is maximal for a pure Gaussian 
state. Further, the standard form (fTS]) implies that for 
every value of | \R\ \ such a state exists, and that the max- 
imal value is given by ||Q||| = 2J2kLi M 2 K| 2 < M/2, 
as |u fc | 2 + \v k \ 2 = 1. □ 

Hence, for every pure Gaussian state with standard 
form (|15p the value of the pairing measure is given by 

M G (p) = 4Efeii M 2 k| 2 - Since \v k \ 2 = 1 - \u k \ 2 the 
measure attains its maximum value for \u k \ 2 = \v k \ 2 = 
1/2, i.e., for the states already identified as maximally 
paired. 

■Mg(p) will appear again when we study the use of 
paired states for metrology applications, linking the pair- 
ing measure to the usefulness of a state for quantum 
phase estimation and giving support to the "resource" 
character of paired states. 



PAIRING OF NUMBER CONSERVING 
STATES 



D. A pairing measure for Gaussian states 



The angular momentum representation of paired states 
depicted in Fig. [1] suggests the introduction of a pair- 
ing measure via a quantity related to | (j k ) PG \ 2 + 

\(j { k V) ) P o\ 2 = \(aW-k) P o\ 2 -- 



3 Maximally entangled states of two qubits share an analogous 
property about entanglement witnesses [4l| . 



In the last section we gave a complete solution to the 
pairing problem for fermionic Gaussian states. There, 
Wick's theorem lead to a reduction of the problem to 
properties of the covariance matrix. For number con- 
serving systems, the situation is more complicated, as 
now also operators of the form ajajofca; have to be taken 
into account. However, we will derive pairing witnesses 
capable of detecting all number conserving BCS states 
in Sec. IV AI using the concept of convex sets. For cer- 
tain classes of BCS states we will construct a family of 
improved witnesses using the analogy to the Gaussian 
states. Witnesses have the drawback that they depend 
on the choice of basis. I.e. even if a witness detects p, 
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FIG. 2: (color online) Expectation values of the vector 
Eq. (|25p . For all number conserving states these lie within 
the convex set CjJ 1 indicated by the dashed back lines. The 
extreme points of the polytope are given by (0, 0, 0), (2, 0, 0), 
(4, 2, 0) and (2, 1, ±1). Unpaired states have expectation val- 
ues in the smaller convex set C^ paircd (solid red) which has 
extreme points (0,0,0), (2,0,0), (4,2,0) and (2, 1/2, ±1/2). 



it does not detect all states related to p by a passive 
transformation. We will show that the eigenvalues of the 
reduced two-particle density matrix can be used to obtain 
a sufficient criterion for pairing in Sec. IV Bl which is basis 
independent. We close the section with the construction 
of a pairing measure in Sec. IV CI 



A. Pairing of all BCS states and geometry of 
paired states 

In a realistic physical setup it may not be practical 
to perform all the measurements needed according to 
Lemma IIII.8I to check the necessary and sufficient con- 
dition for pairing. Having access only to a restricted set 
of measurements, necessary criteria for pairing can be 
derived. In this section we consider the simplest case of 
a symmetric measurement involving four modes, i.e. we 
are looking at the following vector of operators: 



n k + n- k +n t + ri-i 
() : , = [ Tikn- k + niti-i 
a\a}_ k a-iai + h.c 



(25) 



Remarkably, these expectation values will turn out to be 

sufficient to detect all BCS states as paired. 

We are interested in (7" n P airod = {tr(0 3 p) : p separable}, 

the set of all expectation values of O3 which corre- 
spond to separable states. If for some p the vector 
v p = tr(C?3p) is found outside of (7™P aucd then it fol- 
lows from Lemma Hll.81 that p is paired. Membership in 
^unpaired can ^ eag jiy cnec ked by the following Lemma. 



Lemma V.l. A number conserving state p has expecta- 
tion values of O3 (see Eq. i25\) ) compatible with separa- 
bility if and only if tr(H^p) > for k = 1, 2, 3, where 

H i± = 7}( n k + n -k +ni+n-i)- (n k n- k + njra_i) 



± (a\a}_ k a-iai + h.c), 
= (n k n- k + nin-i) ± (aW_ k a-iai + h.c), 



H, 



(p) 

2± 



(26) 
(27) 



H 



3± = I - -{nk + n-k + ni + n-i) + -(nkn-k + nin-i) 



(28) 



are given 



Hence, the extremal points of the set (7™P aired 

by tr(H { k P 2p) = for three of the witnesses ([26 ]) -([28] ) . 
The faces of £7" n P airod consist of points for which at least 

one of the expectation values triH^lp) vanishes. 



h[± and are also pairing witnesses, while H. 

non-negative on all number conserving states. 

Lemma V.2. Every number conserving fermionic state 
fulfills tr(H k± p) > 0, where 

Hi = ~(nk + n-k + nt + n-i) - (nkn-k + nin-t) (29) 
H 2 ± = (n^-k + nin-i) ± (ala'_ k a-iai + h.c.) (30) 

H 3± = 2- -(nk+n-k+ni+n-i)±{a\a}_ k a-iai+h.c.). 

(31) 

The extremal points of the set C^ 1 = {trfO^p) : p E 

O3 

S(A^P) : M,N E N} are given by tr(H k± p) = for 



r(p) 



(p) 
2± 



three of the witnesses (|29p - (|31|) . The faces ofC^ 1 consist 
of points for which at least one of the expectation values 
tr(H k ±p) vanishes. 

The proofs of the two lemmas can be found in Ap- 
pendix [Bl We denote by C ,unpairod and C al1 the poly- 
topes containing all expectation vectors v p correspond- 
ing to unpaired states or all number conserving states, 
respectively. They are bounded by 6 resp. 5 planes de- 
fined through the witnesses given in Lemmas IV. II and 
IV. 21 The situation is depicted in Fig. [2] 

The witnesses h[± given in Eqs. (|26|) allow to detect 
all number conserving BCS states as paired: 

Lemma V.3. The number conserving BCS state 
^scs) given in Eq. |_?7| ) is (except for the trivially un- 
paired cases a k = S kko and N = M) detected by the 
witness Hp by choosing any two modes (k,l). 

Proof. The first two terms in -ffj± are designed such 
that their expectation value vanishes for states such as 

^BCS/- Since we either have a pair or no particles 
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in the modes (k, — k) we are in an eigenstate with 
eigenvalue of the operators n k + ri- k — 2n k n_ k . 
The expectation value of the third term is 
found using the representation Eq. (fT5|) as 

ICjvlWHefafcarJEi^...^-! Kl* ■ ■ ■ \<*m-A 2 , 

which is nonzero unless N — M or all but one ct k 7^ 0. 
The sign can be adjusted by a passive transformation to 



give <flg)g8 fl < 



□ 



This shows that indeed all BCS states are paired, as 
desired. 

The witnesses H^±, while detecting every BCS state 
as paired, are in general far from optimal. As the num- 
ber conserving BCS states appear in many physical set- 
ting, like in the BEC-BCS crossover [42j], it is desirable 
to construction improved witnesses tailored for this class 
of states. For BCS states realized in nature it is often 
appropriate to assume some symmetry of the wave func- 
tion |*scg(afc)) = (Efe=i a fcffc) l°>- For example, if 

- at and if we are dealing 

-fcT k -l 



Rj = at a f - 



P 



kl ' " k+M 

with an isotropic setting, ot k — otk+M will hold. It is 
further often appropriate to assume that the number of 
modes is much bigger than the number of particles, i.e. 
M ^> N. For this kind of states we will construct pairing 
witnesses via the correspondence to the Gaussian picture. 
We sketch the idea of this construction leading to Thm. 
IV.41 and give the details in the appendix. 
We have shown in Sec. Ill Dl the connection of the Gaus- 
sian wave function and the number conserving wave func- 



tion via l^Gauss) = Y,k=i^\^ K BCs( a k))- Consider a 
number conserving observable O and denote by {0)causs 
and (O) n its expectation value for the Gaussian and 2N- 
particle BCS wave function respectively. If the distri- 
bution of I A at 1 2 is sharply peaked around some average 
particle number N with width A, then (O)causs ~ (O) n 
for any integer N € [N - A, N + A] . In Thm. ITVTT1 
we have constructed witnesses H for all Gaussian BCS 
states. As these witnesses are optimal, they suggest 
to constitute an improved witness detecting the corre- 
sponding number conserving BCS state. But H includes 
terms of the form P k that do not conserve the particle 
number. Hence, this witness cannot be applied directly 
to the number conserving case. Using Wick's theorem, 
(PjJ-Pfe+M) Gauss = u k v k {Pl) Gauss holds under our sym- 
metry assumption. This suggests that we replace the 
non-number conserving operator u k v k Pl by the number 
conserving operator P k P k+ M- We define operators 

H k = 2(1 - e - \v k \ 2 )N k - A{PlP k+M + h.c), (32) 
N k = n k + n- k + n k+M + "-(fe+M), (33) 



(TV) 



Theorem V.4. Let M,N G N and let 1 < N < 2M. If 

1 > e > 18/VttN the Hamiltonian H({v k }) — Efc=i Hk 
is a pairing witness detecting 



I* 



(iV) 

BCSySymi 



/ M \ N 

= C N [Y / a k (Pl+Pl +M )\ |0) 



The proof is given in Appendix [Cl 



B. Eigenvalues of the two-particle reduced density 
matrix 

In this section we derive a basis independent condi- 
tion for detecting pairing. The two-particle reduced den- 
sity matrix O contains all two-particle correlations. As 
a change of basis, a\ 1— ► J2 k Ui k a\ leaves the spectrum of 
O unchanged since 

o$ )m - (u® u) (ijUmn) o$ nUpq) (u® u)\ pqUkl) , 

we are lead to the following theorem: 

Theorem V.5. Let p be an unpaired state, and let O 
be its two-particle RDM O. Then \ ma x(0) < 2, where 
\nax denotes the maximal eigenvalue. 

Proof. If p is unpaired, then there exists a separable 
state p a £ S sep having the same two-particle RDM. 
Any separable state is of the form p s = J2 a f-aP^, 
where p^ = \^ a 1)(^\ 7 |^ (q) ) = Ui a U\°) and 
A*a = 1- Here, {a^.}.; denotes some basis of 
mode operators. The RDM is of the form — 
EaMa0 (a) , where is the RDM for the state p( a \ 
The RDM is calculated in the basis {a|}i, and the 
different bases are related by a unitary transformation 
4 = Ej U£>al h , so that 0g> (JU) = tv[p^ a\a]a iak ] = 

(!/ (o) ®P (a) )(«)wOH )W (l/W8^)i, )W , where 



oi a '°\ , n = (at at a„ a„ )„(<»). 



In the basis of 



the {al,.}i the expectation value {a) am a^ an a aq a ap ) p[a) 
is of the simple form (<4 i a a- a a l aa fc )p(<>0 = ^ik^jl ~ 
SuSj k . Hence, the spectrum of the O^ 01 ' is given 
by spec(0( Q )) = {0,2} Va The two-particle RDM 

is hermitian as 0^ ij){kl) = 0(«)(ij) = {a\a\a j a l ) = 
(a\a^aia k ) = 0(ij)(fcj)- Then Weyl's theorem [43| im- 
plies X max (EaMaO") < J2a PaX max {O a ) < J^a 2 A*a < 

2. □ 



where < \v k \ 2 < 1 — e Vfc for e > 0. Further, we intro- 
duce the notation a k = v k /^/l - \v k \ 2 , N = J2 k =i \ v k\ 2 
and we denote by N the biggest integer fulfilling N — N > 
0. Then the following holds: 



An example of a state detected as paired via criterion 
is the BCS-state (T7]) with N = 2, M = 3 and all a k 
equal. The largest eigenvalue of its two-particle RDM is 
given by X max = 8/3. 
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C. A pairing measure for number conserving states 



VI. INTERFEROMETRY 



In Sec. IIVDI we have derived a pairing measure for 
Gaussian states. The correspondence with number con- 
serving BCS states will be a guideline to derive a measure 
for number conserving states. However, the measure of 
Def. IIV.31 involves expectation values of the form (a k a_ k ) 
that vanish for states with fixed particle number. Yet, 
Wick's theorem suggests that a quantity involving ex- 
pectation values of the form (P k Pi) will lead to a pairing 
measure. This is indeed the case, which is the content of 
the following theorem: 

Theorem V.6. Let p be a number conserving pure 
fermionic state. Then the following quantity defines a 
pairing measure: 

M( P ) = 



max < max 



M 



(34) 

where P k — a k a < _ k and the maximum is taken over all 
possible bases of modes {aj}i- For mixed states p, a mea- 



sure can be defined via 

M(p) = min^piMipi), 



(35) 



where the minimum is taken over all possible decomposi- 
tions of p = ^2 i Pi Pi into V ure states pi. 

Proof. The positivity of M and its invariance under pas- 
sive transformations follow directly from the definition. 
It remains to show that M. is zero for separable states. 
We will prove in Lemma ID . 1 1 Appendix [Dl that any sepa- 
rable state of 2N particles fulfills J2ki \( p l p i)\ ^ N > and 
that this bound can always be achieved, which concludes 
the proof. □ 

We close the section by calculating the value of the 
pairing measure for two easy examples. Let 



I* 



BCS I 



k=l 



P 1 



|o>, 



M 



N 



c 



N 



\k=l 



|0>, 



(36) 



(37) 



the tensor product of TV spin-singlet states and the BCS 
state with equal weights, respectively. These states have 
a pairing measure A4(\^ s )) = N resp. M.Q'&gcs) = 
N(M—N). Thus, for the spin singlet the pairing measure 
has in addition the property that it is normalized to 1 and 
additive, while it is subadditive for l^&cs^)- Further, 
this example suggests that the pairing of -MH^bc^) — 
N(M — N) is stronger than for We will see in- 

deed in Subsec. lYl Bl that states of the form \^ e ) allow 
interferometry at the Heisenberg limit. 



The goal of quantum phase estimation is to determine 
an unknown parameter <p of a Hamiltonian H v = pH 
at highest possible accuracy. The value of tp. is inferred 
by measuring an observable O on a known input state 
that has evolved under H v . In a region where the ex- 
pectation value (0(ip)} is bijective, ip can be inferred by 
inverting (0(tp)). In a realistic setup, however, (0(tp)) 
cannot be determined, as this would require an infinite 
number of measurements. Instead, one uses the mean 
value of the measurement results, o, as an estimate of 
(0(p)). This will result in an error Sip for the param- 
eter to be estimated, as for a given value of tp we have 
(0(p)) = o± ^/Var(o). Linearizing around the real value 
of ip, it follows that the uncertainty of p is given by 

HE! 



((Sip) 2 ) = 



Var(O) 
\d(0)/dtp\< 



(38) 



where Var(O) = (O 2 ) — (O) 2 , and we have used the fact 
that Var(O) = Var(o). Further, it can be shown that the 
minimal uncertainty of ip is bounded by [46L l47l | 



((M 2 )Var(tf) > — , 



4i/ ' 



(39) 



where v is the number times the estimation is repeated. 
Eq. ([39|) derives from the Cramer-Rao bound and is 
asymptotically achievable in the limit of large v. 
For a given measurement scheme, i.e. for a given input 
state and a given observable O, the uncertainty in tp can 
be reduced by using N identical input states and average 
over the iV measurement outcomes. As the preparation 
of a quantum state is costly, a precision gain which has 
a strong dependence on N is highly desirable. If these 
probe states are independent of each other, the preci- 
sion scales like 1/V~N. This is the so-called standard 
quantum limit (SQL). Using distinguishable or bosonic 
systems, this limit can be beaten by a factor of \HV 
by using number-squeezed input states [H, HH , 

iV-particle NOON states or maximally entangled GHZ- 
states -^(|iV,0) + |0, N)) |H H HE H3- Achieving 
this so-called Heisenberg-limit is the big goal of quan- 
tum metrology. 

Less is known for fermionic states where number squeez- 
ing and coherent TV-particle states are prohibited by 
statistics. Nevertheless, there exist fermionic iV-particle 
state which can achieve the Heisenberg limit for phase 
measurements in a Mach-Zehnder interferometer setup 
[5^ . Taking the existence of such states as a start- 
ing point, we show that paired fermionic states can be 
used as a resource for phase estimation beyond the SQL. 
We will consider two different settings. The first set- 
ting will be the standard Ramsey-interferometer setup 
of metrology, where the coupling Hamiltonian is propor- 
tional to the number operator. Here, we will see that 
paired states lead to a precision gain of a factor of 2 
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I in/ 




a, ,a_ , 



In the next section we will derive the best possible pre- 
cision obtainable by using unpaired states, and compare 
this result to the precision achievable by using paired 
states. It will turn out that already at two-particle level 
paired states have more power than the unpaired states 
for our setup. 



FIG. 3: Scheme of the Ramsey interferometer setup. The 
incoming wave function \^i„) enters a beam splitter (BS). 
Then particles in the modes a) ±l evolve under the Hamiltonian 
ipH . At the end a particle number measurement is performed 
on all particles. 



compared to separable states. The second setup involves 
a more complex coupling. Here it will turn out that by 
using paired states the Heisenberg limit, i.e. a phase sen- 
sitivity (Sip) 2 ~ 1/N 2 , can be achieved. 



A. Ramsey interferometry with fermions 

1. General setup 



We consider the standard Ramsey interferometer setup 
where a state in the modes {a\.., c4-}jf=-M undergoes 
mode mixing at a beam splitter, 



H 1 I 



a 



't 

±U 



(40) 
(41) 



before evolving under the action of the Hamiltonian 

M 

H N =J2(n l} +n. h ). (42) 
i=i 

Finally, a particle number measurement is performed on 
the system, to compute the parity 



■p = (_X)Z^j "0 +"l 



(43) 



where - "' ""' "' ' 



n o = a kj a kj an d n i' = a '- kj a '~ kj - According to 
eq. (|38p . the phase sensitivity is given by 



{6<pf 



1 - (Vf 



dip 



(p) 



(44) 



where we have exploited V 2 = 1. Due to the fermionic 
statistics the parity operator can be written in the form 



M 

V = n (l - 2(r# ) + nf) + 4n^n^) 



(45) 



2. Bound on unpaired states for the standard interferometer 

In this section we derive a lower bound on the phase 
sensitivity when using an unpaired state of 2N particles 
as input states. 

Theorem VI. 1. For the Ramsey interferometer de- 
scribed above the phase sensitivity is bounded by 



(^) 2 > 



1 



2vN' 



(46) 



when an unpaired state of 2N particles is used as input 
state. 



Proof. We will use (|39|) to derive the bound. Hence, 
we have to estimate an upper bound for the variance 
of the Hamiltonian Hm dchncd in (|4"2")l . As Hn as 
well as H^j contain operators from the set A 2 only, it 
is sufficient to proof the bound for product states, as 
for every unpaired state there exists a product state 
having the same expectations. In Lemma IA.2I of Ap- 
pendix [X] we have shown that for pure separable states 
{n k n t ) = \P kl \ 2 - P kk P u + P kk 5 kh where P € C 4Mx4M 
is a projector of rank 2N. We arrange the indices as 
—Imi ■ ■ ■ i hii ~kM, ■ ■ ■ , fcjvr and partition the projector P 

Midi that P = I f+ ^ I , where A,B,C G € 2Mx2M . 



fit c 



Then 



2M 

Var(H N ) = Yl A > 



2M 

£ 



-4; 



tr^BB 1 



(47) 



as Hn only involves the modes —Im---^m- In the last 
step we have used P 2 = P , implying A — A 2 — BB^ . As 
rank(P) = 2N : there exists some unitary U such that 

12N 



P = UldwW, where Id 



2N 







G C 



4M x 4M 



Partitioning the unitary [/ = I l^ 11 H 12 ) , where Ua G 

L'21 <^22 



C 



2Mx2M 



J = l 



«, j = 1, 2, the projector P is of the form 



U21M2NUI U 21 ld 2N U^ 



Using the the above representation of P and the cyclic- 
ity of the of the trace, we can write Var(i7/v) = trL4£?] 
with hermitian matrices A = Id2Ar/7-J 1 /7nId2Ar, B = 
Id2ArL^2it^2iId2Ar. The trace can be interpreted as a scalar 
product maximized for linearly dependent A and B. Ex- 
ploiting the unitarity of U, one sees immediately that 
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the variance is maximized for A — c/(l + c)Id27v, B = 
1/(1 + c)Id2Ar for some constant c. Hence, Var(-ffjv) < 
c/(l + c) 2 tr[l 2 7v] < N/2. Inserting this into eq. ([39]>. we 
find that (Sip) 2 > □ 



3. Inter ferometry with two particles 

In this section we will show that already a two-particle 
paired state can beat the bound for the phase sensitivity 
using unpaired states (|46| . Hence pairing manifests itself 
as useful quantum correlation already at the two-particle 
level. We show the following: 

Theorem VI. 2. Using the paired state 

M 

l*& } > = ( E a i<aU I |0) ' (48) 

with normalization Y] ^ \cx.j\ - 1 " s input state for 

the Ramsey interferometer, the optimal phase sensitivity 
is given by 



(Sep) 



2 

min 



2(l + 2E" 1 ReM J ) 



1 

> -. 



(49) 



Proof. Take \^ in ) as the input state. After an appli- 
cation of the beam splitter transformation (j40|) and an 
evolution under the Hamiltonian (|42|) . the measurement 
outcome of the parity operator is calculated to be 



(V) = 1 - sin 2 V\1 + 2J2 Rc(oyft 
i=i 



(50) 



Using eq. (|44|) . we obtain (j49|) . The bound of 1/4 can be 
obtained for a state where a% — (3k V k. □ 



Thm. IVI.2I shows that there exist two-particle paired 
states exceeding the bound on product states PS)) . 



4- Inter ferometry with 27V -particle BCS states 

Generalizing the result obtained in the last sec- 
tion, it follows immediately that states of the form 
\^m)® N i wm l ea d to a phase sensitivity (Stp) 2 ^ = 



In this section we will 



1/ [2JV(l + 2ESiBfi(a J 4-)) 
show that the same result can be achieved using BCS 
states. 



Theorem VI. 3. Let the paired state 

I M 



where we use the normalization condition J2j\ a j\ 2 + 
\f3j\ 2 — 1 be the input state for the Ramsey type interfer- 
ometer defined above. Then the optimal phase sensitivity 
is given by 



(Sip) 2 



1 



2AT(l + 2£,Re( ai /3,-)) 



(52) 



Proof. Like in previous sections we will use the corre- 
spondence to the Gaussian state 



I* 



(2ft) v 
in. Gauss I 



cexp 



M 

J=l 



10), 



where — jV| <C iV for the calculation. After the 
state has passed through the interferometer, the expec- 
tation value of the parity operator is readily computed 
to be (T) G auss = Uji 1 ~ \c?\a 3 + f3 3 \ 2 sin 2 p). As only 
number operators are involved, (V) Gauss ~ ('P)n, where 
(. . .)jv denotes the expectation value of V for the state 



)• Expanding (|4"4")l for small value of tp and us- 
ing N = \c\ 2 J2k l«fc| 2 + \Pk\ 2 = \c\ 2 , one obtains (52]), 
which has minimal value (Sip) 2 = l/(4iV). This result is 
obtained when ctj = f3j V j. □ 

The above result shows that paired states result in a 
precision gain of up to a factor of 2 compared to the best 
precision obtainable for unpaired states (|4"5)) . 
The pairing measure derived in Sec. HVDI and lVCI auan- 
tifies the precision gain obtainable by the use of paired 



states. To see this, denote by |* 



(2N)f 
in, Gauss l 



the state af- 



ter the beam splitter transformation. Then the pairing 
measure pV.3|) for this state evaluates to 



Mg{\<1% uss )) 



N 2 



l + 2^Re(a J /i i ) , 



so that 



(M 2 



N 



(53) 



The above relation demonstrates that M. is indeed 
quantifying a useful resource present in paired states. 
Whether this interpretation can be extended to mixed 
states will not be explored here. 



B. Interferometry involving a pair-interaction 
Hamiltonian 



So far we have seen that paired states lead to a gain 
of a factor of 2 in precision compared to unpaired states 
in a Ramsey-type interferometer. This section will show 
that paired states are even more powerful and can lead 
to a precision gain of a factor of N when measuring the 
phase of a pair-interaction Hamiltonian. 
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bl 



iH<p 



PmPi 



FIG. 4: Setup which allows interferometry with paired states 
at the Heisenberg limit. Particles in modes a\. and b\. evolve 
under the complex coupling Hamiltonian H (for the detailed 
form of H refer to the text). In the end particle numbers are 
measured. 



We consider a setup where two fermionic states enter the 
ports A and B of an interferometer. The particles enter- 
ing port A can occupy the modes {a\} k L_ M , while the 
particles entering through port B can occupy the modes 
{bt}kL-M- Then the two states evolve under the Hamil- 
tonian H c to be defined below and a particle number 
measurement is performed at the end. The situation is 
depicted in Fig 2J We will compare the power of paired 
states over unpaired ones for two different settings. Wc 
start by introducing some basic notation: 



Prerequisites 



We define pair operators Pi = a\a}_ k and Q\ = b\b*_ k 
and their equally weighted superpositions 



M 



M 



Pm 



<M 



£4 



fc=i 



'M 



(54) 



with normalization constant cjy = {NM\/M N ) 
while the separable states are given by 



1$ 



\^ 2N) )a\^ N) )l 



(59) 



where \4>^ 2N ^)a,b are separable states in the modes a k and 
b k respectively. 

After the input state has evolved under the Hamilto- 



nian H v into the state 



observable O is used as an estimator to determine the 
parameter ip to a precision given by ([55)) . Instead of 
working in the Schrodinger picture of state evolution it 
turns out to be more convenient to tackle the problem 
in the Heisenberg picture, where O evolves according to 
O —> O' = e -iHa V Q e iH cV _ We are interested in the phase 
sensitivity for small ip, so that we can expand (|38|) in 
powers of ip, arriving at 



0(<p) 



= 0-i<p[H c ,0] 



1 



-<p 2 {H 2 + OH 2 -2H c OH c ) 

(60) 



If the input state l^jy ) is an eigenvector of O with 
eigenvalue 0, we obtain the following simple expressions 
for (O) and Var(O) 



^-(O) = A^\(H c OH c )\ 2 +0{^), (61) 
Var(O) = ^ 2 (iJ c 2 iJ c ) +C(^ 3 ), (62) 
so that the phase fluctuation (Sip) 2 simplifies to 



fe=i 



The operators p M and q M fulfill the commutation rela- 
tions 



1 



Pm>Pm] = -1 + -rjNa 



Wm^m] = -1 



M 

1 

w 



-N h 



(55) 
(56) 



where rik = a k a,k so that N a = J2k( n k + n -fc), and 
N b = J2k( m k + m -k) with n k = a\a k and m k = b\b k 
being the number operators for particles in modes a k and 
b\ respectively. 

We will compare the power of two paired states and 
two unpaired states entering through port A and B. The 
bound for unpaired states will be derived again via (f39|) . 
As H c and H 2 will be elements of Ai , it is like in the last 
section sufficient to compare the power of paired states 
to those of sparable states. The paired states will be of 
the form 



ffW) = \N)[ M W b M) , 



■ N 

iao1 m) = 4 m Vm)», iA0i M) = 4 A %ino>, 



(57) 



(58) 



(M 2 



(HQ 2 H) 
4\(HOH)\< 



0{!p) 



(63) 



An observable fulfilling this property is O = [n]^) 2 , 
where n { ~ ] = \(p* M p M ~ OmQai)- 

The commutation relations for p\ d and gL (|55p imply 
that in the limit of infinitely many modes M — > oo the 
operators p M and (f M become bosonic. We will thus start 
out with a scenario where the input states are in the 
bosonic limit and then turn our attention to a setting 
which is far from the bosonic limit. 



2. Bosonic limit 

In this section we will consider the scenario M — > oo, 
i.e. we are in the bosonic limit, where the limit is taken 
for the expectation values of the operators. We will con- 
sider a coupling of the form H c = ipHryo where 



(64) 
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and measure (n^) 2 . 

We start deriving the best precision for unpaired states 
using (j59"j) . We will use a finite M for input state, 
coupling Hamiltonian and measurement and then take 
the limit M — > oo. To be precise, the calculation will 
be done for H M — ^(a\ 4 b M + om&m) and (n^) 2 . 
Then limM^oo(^a,6|-ffiw|^a,b) = due to the conserva- 
tion of particle number. Hence, limM— »oo Var(_ff^/) = 

]imM-»oo(-ffAf) = ^m^oc((PmPm) + (Qm9m)) 2 = 0, as 

(pIpm) = jjEki \( p l p i)\ 2 < N / M i where the last in - 

equality results from the bound of the pairing measure on 
unpaired states Thm. IV. 61 The same holds for (q M qM)- 
Hence, in the limit M — > oo the variance of vanishes. 
For the setting of paired states, however, we can obtain 
the following result: 

Theorem VI. 4. For paired input states, the interfer- 
ometer depicted in Fig [7] allows to estimate the coupling 
parameter tp to a precision 



i 



2N 



2 ' 



(65) 



Proof. Consider an interferometric setup depicted in Fig. 
IU where the 2iV-particle input state and the coupling 
Hamiltonian are defined in eqs. (|57|) and (|64|) respec- 
tively. We will again use a finite M for input state, 
coupling Hamiltonian and measurement and then take 
the limit M -> oo, i.e. we use |* in > = \N)\^ ] \N)[ M \ 

Hm = Kaj^&M + «m&m) an d {n^) 2 . Making use of 
the relations 



p M \N)W = a N \N-l)i M ^ 



(66) 
(67) 



where ajy = \/N(l - (AT - Tj/M), a lengthy but 
straightforward calculation leads to (8tp) 2 M = a 1 — 3- + 

0(tp) using (j6"5j) . Taking the limit M — > 00 leads to the 
result of the theorem. □ 



3. Interferometry far from the bosonic limit 

In the preceding section we have studied the power 
of paired states in the bosonic limit. As the power of 
bosonic particles for interferometry has been known for 
quite a while, the use of paired states where the fermionic 
nature of the particles survives might be a more interest- 
ing question. In this section we will show that even far 
from the bosonic limit paired states can achieve a preci- 
sion gain of order N for quantum metrology. 
We will study a coupling Hamiltonian of the form H c = 
tpHp where 



H F =ypiQ k +p k Qi 



(68) 



First, we will give a bound for the phase sensitivity 
achievable by using product states at the input: 



Theorem VI. 5. Using product states of 2N particles 
as input states for the interferometric setting depicted in 
Fig. the phase if of the coupling Hamiltonian H c = 
Hp, where Hp is defined in (|68[) can be measured to a 
precision (5ip) 2 > l/(16Af). 

Proof. For every product state of the form (Hp) = 
due to particle number conservation. Hence, Vai(Hp) = 
(Hp). We will bound this expectation value: 

(Hi) = Y,{Pl p i){Q\Qk)+c.c. + Y,(PlPk)(QkQl) 

kyLl fc 

I v 1/2 1 \ 1/2 

k 

/ . xl/2 

From Lemma ID. ll we know that {J2k^i\( P k Pl )\ 2 j — 
VN. Further, (ffcP fc t) = (1 - (n k - n_ k ) 2 - n k n_ k ) < 1 

and J2k( P k P k) < N - Thus Var(iJ p9 ) < 2y^Ny r N+2N = 
AN which leads immediately to our result via (|39|) . □ 



This bound can be beaten by a factor of vA" using 
paired states. A lengthy but straightforward calculation 
leads to the following result: 

Theorem VI. 6. Using paired states of the form (|57j) 
as input states for the interferometric setting depicted in 
fig. [7J the phase tp of the coupling Hamiltonian H c = 
<pHp, where Hp is defined in f|68|) can be measured to a 
precision 



(Stpf 



M(M - 1) 



8N(M - N)(M - 1 + MN — N 2 ) ' 



(69) 



This theorem implies (Sip) 2 ~ 1/N 2 for all M > 2N. 
In conclusion we have shown that paired states are a re- 
source for quantum metrology. Theorem I VI. 6 1 is the main 
result of this section. We have remarked already at the 
beginning of this chapter that it has been proven before 
that the Heisenberg limit can be achieved using fermionic 
particles [55(. However, these states were constructed in 
an abstract way, while we prove that the BCS states that 
can be created easily in an experimental setup are a very 
powerful resource for quantum metrology. 



VII. APPLICATION TO EXPERIMENTS AND 
CONCLUSION 

In summary, we have developed a pairing theory for 
fermionic states. We have given a precise definition of 
pairing based on a minimal list of natural requirements. 
We have seen that pairing is not equivalent to entan- 
glement of the whole state nor of its two-particle re- 
duced density operator but represents a different kind of 
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quantum correlation. Within the framework of fermionic 
Gaussian states we could solve the pairing problem com- 
pletely. For number conserving states we have given suf- 
ficient conditions for the detection of pairing that can 
be verified by current experimental techniques, e.g. via 
spatial noise correlations [H, H3, [si] and prescribed a 
systematic way to construct complete families of pairing 
witnesses. 

To shed some light on the pairing debate 0, E3, EH, 
[IH, we would need access to the proportionality fac- 
tor linking the quantity plotted in Fig. 4 of 0] to the 
local pair correlation correlation function G2 {r, r) — 

(*{(r)*{(r)tt T (»•)*!(»•))._ 

Another important point of our work is the utility 
of fermionic states for quantum metrology. While it 
has been shown that, in principle, fermionic states can 
achieve the Heisenberg limit for precision measurements 
in a Ramsey-type interferometer 55| , we could prove the 
usefulness of states that are available in the laboratory. 
Furthermore, the optimal precision for the Ramsey-type 
setup is proportional to the pairing measure introduced 
from an intuitive picture in Sees. Ill CI and |Vj This en- 
dows the measure with an operational meaning. 
The results we have presented are just a first step in 
understanding pairing and its relation to other types of 
quantum correlations. 

We hope that the pairing theory we have developed will 
help to get a better understanding of correlated many- 
body systems, and can provide a new perspective on 
quantum correlationsa and may serve as a starting point 
for further inquiries. 

For example, one might attempt a finer characteriza- 
tion of pairing, e.g., Y%=i P k\°) and Sfeii P kW repre- 
sent paired states of rather different nature: it would be 
interesting to develop witnesses or measures which allow 
to determine over how many modes the pairs in a given 
states extend and to relate these differences to applica- 
tions in metrology ore elsewhere?. Moreover, the the- 
ory we developed has been concerned with finitely many 
modes only and it is an obvious question whether gen- 
eralizing to an infinite dimensional single-particle space 
gives rise to new phenomena. 

Up to now we have concentrated on fermionic states. 
But the question of pairing in bosonic systems might be 
equally interesting and relevant for recent experiments 

What about higher-order correlations? The set of un- 
paired states contains both separable and highly corre- 
lated states. This is, for example, reflected in the fact 
that there are unpaired states which can be transformed 
to paired ones by single-mode particle number measure- 
ments (e.g., (a\a\al + a\a\ag)\0) by measuring particle 
number in mode b — 0,3 + a&). A theory of higher-order 
correlated states could be developed along the lines dis- 
cussed here, e.g., by changing the set of observables on 
which the states are compared to uncorrelated ones and 
defining as nth-order correlated those states whose ex- 
pectation values on nth-order observables cannot be re- 



produced by 771 < n-correlated states. 

Tools and methods from entanglement theory have 
been very useful in analyzing pairing. One very impor- 
tant such tool, however, is missing: positive maps, that 
is, transformations which do not correspond to a physical 
operations but nevertheless, when applied to a subsystem 
in a separable state with the rest, map density operators 
to (unnormalized) density operators and thus provides 
strong necessary conditions for separability. Finding an 
analogy might prove very useful for the analysis of many- 
body correlations. Another importat object in the the- 
ory of entanglement is the set of LOCC operations (local 
operations and classical communication), i.e., the oper- 
ations that cannot create entanglement. In the case of 
pairing, the analogous set would contain passive oper- 
ations and discarding modes. Are there other physical 
transformations that cannot create pairing? Do paired 
states, then, possibly allow to implement such transfor- 
mations similar to entanglement enabling non-LOCC op- 
erations? 
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APPENDIX A: USEFUL PROPERTIES OF 
SEPARABLE STATES 

We give two technical lemmas which involve useful 
properties of product states. 



1. Bound of a]a^akai + h.c. on separable states 

In this section we prove a bound of a some special two- 
body operator on product states: 

Lemma A.l. Let p G S sep be a separable state. Then 



\tv[{a\a\a k ai + h.c.)p s }\ < 1/2 



(Al) 



Proof. Let Kya be the Hilbert space spanned by 
a\,aj,al.,aj and define Ayjy = ajaj-afcdj + h.c. 
Then tr[A ijkl p] = tr[p^A ijkl ], where p^ k » = 

En=o Pijkl\ n )( n \ * s a rnixed separable state accord- 
ing to Lemma III. 21 and \n) denotes the occupation 
number basis for the subspace Tiiju- It is eas- 
ily checked that Aijki can have non-vanishing expec- 
tation value only for the two-particle state |2) 



,k,l ^r 11 



(V. |0). Using |2Rc(a6)| < 
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\a\ 2 + \b\ 2 for any complex numbers a, b and the nor- 
malization conditions X)r l/^l 2 = Srl^l 2 = •"■> one ar ~ 
rives at |trL4y fc i/o]| = 2|Re[(/Xii/j — fJ,ji/i)(ukVi — MVk)* 



(iMii a +i^r)(i^i a +iwi a )+(kii 2 +kji a )(i^i 2 +hi a ) 

0.25 + 0.25 = 0.5. 



< 
□ 



2. Expectation values of one-and two-body 
operators for separable states 

In this section we will prove that the one-and two-body 
operators for separable states can be expressed in terms 
of matrix elements of projectors. 

Lemma A. 2. Let p £ S^J be a pure separable state. 
Then 



(rii) = P u 

{a\a\a k ai) p = (P <g> P)(y)(; fe ) - (P <8> P) («•)(*!), 



w/iere P = P 2 = pt is a projector of rank N . 



(A2) 
(A3) 



Proof. Consider M modes. We go into the basis where 
the pure separable state is of the form |$) = Ylf =1 a^jO). 
In this basis (a^a^o^ai)!*) = $itfjk ~ $ikSjh i- e - EH 
for P = Idjv, where Idjv = 1at0O M -jv € <D MxM . Now 
let aj = E fe ^4 fe - Then 

(a|atojfeOj> p (jf) = (U Id Ar J7 t ) ® (C/Id A r?7 t ) ( y )(;fc) - 

= (P®P)(y)(ifc) ~ (P® P) 
and P is a projector of rank AT. 

For the one-particle operators, we obtain (rii) = Pa, as 
(N - l)(m) = ^(rnnj) = Y, jfi Pi/'ii - I Pi 1 2 = 
V . PiP^ - /',.,/',„ = tr[P]P l4 - {P 2 ) u = {N- 1)P U . □ 



APPENDIX B: PROOF OF LEMMAS IVTT1 AND 

El 

1. Proof of Lemma IV. II 

Proof. As , P^± i are built up of operators that 
are the product of at most two creation and annihilation 
operators, we can prove the lemma for separable states. 
In the first step, we will show that the three operators are 
positive on all separable states. Then we will show that 
all states within the set bounded by Ll[± , H%± , Pg± cor- 
respond to a separable state. Finally we will show there 
exist states that are detected as paired by h[± and Hg± . 

Positivity of on all number conserving states will be 

shown in the proof of Lemma IV. II following below. 

To show positivity of h[ ± , , P^z± it is sufficient to 



show the positivity for pure separable states, as the re- 
sult for mixed states follows from convexity. From now 
on, let p e Ssep ■ 
trpgd > 

In Lemma IA.2I (see appendix [AJ we have shown that the 
expectation values of number conserving one-and two- 
body operators can be expressed in terms of matrix el- 
ements of projectors. Let P be the rank N projector 
such that (afata*aj) p = (P <g> P) - (P <8> P)(ij)(kl) , 
(m) = Pa and let P = P\k-k,l,-l the 4x4 principal sub- 
matrix of P where the indices run over k,—k,l,—l. Then 
we have the following inequalities: 

(n k n- k ) = P kk P- k - k - |Pfc_fc| 2 

< \{\P kk \ 2 + \P- k - k \ 2 ) ~\P k - k \\ (Bl) 



\{a\a ] _ k a-iai + h.c.)\ 

= 2|Re(PyP_ fe _ ; - P k -lP- H )\ 
< 2(|p y ||P_ fe _ ; | + |p t _ i ||P_ w |) 



< (IP 



kl 



iP-k-lf + lPk-if + lP-kll 2 ) 



These results imply 



tr[pP^]>-tr[P-P 2 ] 



1 



|P fc - fc | 2 



\Pl-l\ 2 



(B2) 
(B3) 

(B4) 



We use the inclusion principle [43j , stating that the eigen- 
values of a r x r principal submatrix M r ofanxn Her- 
mitian matrix M fulfill A fc (M) < \ k (M r ) < A fe+n _ r (M), 
where the eigenvalues are arranged in increasing order. 
As P is a projector, we have < X k (P) < Afc(P) < 
Afc +M -r(P) < 1- Hence, 

tv[H<gp] > itr[P-P 2 ] > ~5>*(P)(1- A fe (P)) > 

k 

(B5) 



l±a\a}_ k a-iai + 



trp^p] > 

Define 0\ = n k n- k + mn—i > 0, 
h.c. > 0. Then H$ = O x O% , and as [O u Of] = we 
conclude that H$ = O x O^ > 0. 
tv[U^p] > 

We will need the lemma HU] tr[H$p] = tr[H^2p k i}, 

where p kl = J2t=o Pn\n){n\, (3 n > 0, ^ =0 A = 1 and 
\n),n — 0,...,4 are separable n-particle states. Let 



(H 



(p)\ _ 

3±/« — 



(n\H. 



(p) 

3± 



n 



Then a straightforward calcula- 



tion leads to (H^)q = 1 
r(p)\. - n „„a /w(p) 



(H. 



(p)\ _ l 



3±/l " 2' \"3±/a — 2' 

H^)s — and {H^}^ = 0. Linearity of the trace im- 



plies tr[H^p] > 0. 
Hence, all separable states lie within the set bounded 
by the planes defined by the witness operators 

tt(p) tt(p) tt(p) 
i3 l±'- n 2±' i7 3±- 

Next, we show that each point within the polytope 
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^unpaired correS p 0nc j s to a separable state. As S sep is con- 
vex, it is sufficient to check that for every extreme point 
of (7 un P aired there exists a separable state. This is indeed 
the case: The extreme points of C unpairod are (0,0,0), 
(2,0,0), (4,2,0), (2, 1/2, ±1/2) which correspond for ex- 
ample to the separable states |0), a k a\ |0), a k a_ k a l a_ l |0) 
and (at + a |)( a i^ ± al z )/2|0) respectively. 
It remains to show that and are pairing 

witnesses. Define |*) = -±=(a\a^ k + ajal ; )|0). Then 

tr[^|*)<*|]=tr[^l*>(*|] = -l. □ 



2. Proof of Lemma [VT21 

Proof. It is sufficient to prove the lemma for p S S(An), 
as the result for a general number conserving state fol- 
lows from convexity. 
tr[Hip] > 

Hi = \{n k - n_ fc ) 2 + \{ni - n_i) > 
tr[H 2± p] > 

Shown in the proof of Thm. IV. 21 
tr[H 3± p] > 

Let p k i — J2n=l Pn\0)(n\ be the reduced density opera- 
tor in the modes ±fc, ±L We can rewrite H 3 ± in the form 
H 3 = 2-\(n k +n_ k + n l +n_ l ){l^(a\a}_ k a_ l a l + h.c.)). 
Defining C^zp = 1 =F ( a \ a ^- k a -i a i + h.c), we obtain 
(0 2T )o = (0 2T )i = (0 2T ) 3 - (0 2T ) 4 - 1,(0 2T ) 2 < 2. 
This implies tv[pH 3 ±] > 4 - + 2/3 2 + 3/3 3 + 4/3 4 ) > 

4 - H«=o n Pn = 4 - tr[(^fc + n k +n t + n-i)p k {\ > 
As in the proof of Lemma IV. II it remains to show that 
the extreme points of C al1 correspond to some fermionic 
state. It has been shown in the proof of Lemma IV. II 
that (0,0,0), (2,0,0) and (4,2,0) can be reached by 
some separable state. The remaining two extreme points, 



Due to the particle number constraint ^2 k Li(^k) = N 
this value can be taken for k — l,...,N/4. Hence, 



(2, 1, ±1) correspond for example to the state (a k 



)|0). 



□ 



APPENDIX C: PROOF OF THM. IVTil 

In this section we will provide all the details leading to 
Thm. IV. 4[ starting with the bound on separable states: 

Lemma C.l. Let p € S sep with tr[N op p] = N and let 
H({v k }) be as in Thm. \ V.4\ Then tr[H({v k })p s ] > 0. 

Proof. The operator H k acts non-trivially only on the 
modes a\, a)_ k , a\ +M , at_, k+M y Denote by p k the re- 
duced density operator obtained when tracing out all 
but these four modes. According to Lemma [11.21 p k = 

^2n=o Pk l n ) ( n l i s a convex combination of separable n- 
particle states |n)(n|. We proved in Lemma I A. 1 1 that 
Itr^Pfe+M + h.c.)p s ]\ < 1/2. Hence (H({v k })) > 
2N(l~e)~2j:ll 1 (\v k \ 2 (N k )+p k 2} ). Now2/3 fe < (N k ) < 
4 and \v k \ 2 < 1 - e so that |^ fc | 2 <iV fc > + [3 { k ] ) < 4. 



-2Efcli(kfc| 2 (^fc) + 0£>) > -2 • 4(1 - e)N/4 so that 
(H({v k })) > 0. □ 

To show the witness character of H({vk}) we also have 
to prove that there exists a BCS state that is detected 
by the Hamiltonian. We will need the following theorem 
about the distribution described by the | A^v | 2 in (fT6| : 



Theorem C.2. Let \^ Gauss ) = Ew=o A 



N 



BCSI 



like 



m, (53- IfEkLi Wk\ 2 \v k \ 2 = O(N^) for some 7 > 

oo the \Xn\ 2 converge to a 



0, then in the limit N 
normal distribution, 



i 



27TCT a 



exp 



JV 



(N - Nf 
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(CI) 



where 2N = 2j^ fc=1 \v k \ 2 is the mean particle number 
for the variational state, and the variance is given by 



N 



4Ef=ikl 2 KI 2 - 



Proof. For the proof we will need a theorem from proba- 
bility theory known as Lyapunov's central limit theorem 
0: 

Theorem C.3. (Lyapunov's central limit theorem) 
Let Xi , X2 , . . . be independent random variables with 
distribution functions F\,Fi, ■ ■ ■ , respectively, such that 
EX n = p n and VarX„ — a 2 t < 00, with at least one 
a„ > 0. Let S n = Xi + ■ ■ ■ + X n and s n — ^/Var(5„) = 
\J o~\ + • • • + (J 2 ■ Lf the Lyapunov condition 



2+8 

Sn k=i 



E\X 



k — Mfc 



| 2 +<5 



is satisfied for some 5 > then the normalized partial 
sums l5 '"~" E ' S " converge in distribution to a random vari- 
able with normal distribution iV(0, 1) 

Consider the observables X k = n k +n- k , k = 1, . . . , M, 
where n± k = 0, 1 is the number of particles with quan- 
tum numbers ±fc respectively. The X k can be considered 
as classical random variables since they commute mutu- 
ally. In the variational BCS-state the random variable 
Sm = Ylk=i -^-k is distributed according to the probabil- 
ity distribution 



P(S M = 2A0 



A I 



E l Ufc i I 2 ■ ■ ■ | 2 K N+1 1 2 ■ ■ ■ Km I 2 



fcl<...<fcM- 1 



n 



Uk\ 



M 

E 



k / k±<...<kM 

\c\ 2 



=1 \u kl \ 2 ■ ■ -\u kN [< 



{N\) 2 \C N \ 



I A 



(AO 1 2 
N I ■ 



(C2) 



(C3) 
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With the help of Thru. IC.3I applied to the random 
variable Sm we can now complete the proof of Thm. 

IC.2| i.e. show that X^ 1 ^ converges to a normal distri- 
bution for large M. We start calculating the expectation 
value /it of X k . For a BCS-state, X k = 0, 2, as particles 
with quantum numbers ±fc always appear in pairs. As 
P{X k = 0) - \u k \ 2 ,P(X k = 2) = M 2 , we get fi k = 
2\v k \ 2 and E(S M ) = 2£\. K| 2 = 27V. For calculating 
the variance, note that X k = n k + n 2 _ k + 2n k n_ k — 0, 4, 
and P{Xl = 0) = |u fc | 2 , P{Xl = 4) = \v k \ 2 . Hence 



where 



Var(X fe ) = 4K| 2 H 2 , 



4 f = 4j>*|V| 2 . 



A- 



(C4) 



To apply the central limit theorem, we consider 
E(\X k - ^| 4 ). Using P(\X k ~ ^l 4 = 4) = |u fc | 2 , 
P(\X k - fi k \ 4 = (2 - Mfe ) 4 ) = M 2 , and /i fc = 2| Ufc | 2 
we arrive at 



A I 



H = -2j2\v k \ 2 N kl 

k=l 

w = 2j2HPk+M+pt +M p k 



We start with a bound for \{H ) var - (H ) N \ < T x + T 2 , 
where 



Ti 



E iwi^^o) 7V+A - (Hq)n) 
E (IAjv+aI^^jv+a - (^o)iv; 



(C7) 



A bound for T 2 can be easily derived noting that 



E{\X k -ix k \ A ) = 16K|> fc | 2 (|u fc | 6 + |w fc | 6 ) 
< l&\u k \ 2 \v k \ 2 (\u k \ 2 + \v k \ 2 ) 
= 16\u k \ 2 \v k \ 2 . (C5) 

Setting S — 2 in the Lyapunov condition, we obtain 



(Ho)n — (Ho) n ' \ — 



E \ v k\ 2 {{n k ) n - (n k )n>) 



(C8) 



and for n = A + A >= A we have ^2 k |ufc| 2 ((n.fc)„ 
(rifc)jv) < Lfe M 2 ("fc)™ < as |w fc | 2 < 1. Hence, 



M 



1 

3. 



^£;(|A fc - Mfe | 4 ) 



< 



0(JV" 



0, (C6) 



M fe = 



fc=l 



where we have applied the assumption of the theorem 
Efcli KI 2 KI 2 = 0(A 7 ) in the last step. The central 
limit theorem implies that Sm converges to a normal 
distribution with expectation values 2A = 2 J2 k \ v k\ 2 and 
4£ fe M 2 K| 2 . " □ 



variance er 



With this result at hand we can prove the following: 

Lemma C.4. Let H({v k }) and \^BCs) ^ e defined as in 
Thm. [73 Ife > 18/VttN, then 

(y { BCS,syJ H ({ v k})\y { BC S , sy J < 0. 

Proof. We will use the correspondence of variational and 
number conserving BCS states, deriving first a bound 
for \(H({v k })) var - (H({v k })) N \, where |*g „ ss ) = 

El^^UsyJ with \<k sy J ^ke in Thm. EH 
To do so, we will need that the |A„| 2 are normally 



distributed. From \u k \ = 1 



I 2 > e where e > 



18/(V2ttA) and Y?k=i M 2 = N, it follows that a N = 



J2k=i \ v k\ 2 \uk\ 2 > eA = 0(\/A). Hence, we know from 
Thm. IC.2I that the |A n | 2 describe a normal distribution 
around A « A with standard deviation er^y . 



Now, write H({v k }) = H — 2W + 2(1 - e) J2 k =i N k, 



T 2 < 16 



E |A.v + a| 2 |A + A| 



< 8^6-^/2 



2tt 



4A(1 - Erf(cr 



A? 



<8-^+4A(l-Erf(^/y2)), 



(C9) 



where we have approximated the sum by an integral in 
the second step. For bounding Ti, we will show first that 
for n — A + A where A S [— (7 %,o' 2 ^] we have (n k ) n — 
{n k ) n -i > 0. Expanding the BCS wave function 



2M 



< l " ) ^ a 31 ...«,„4...p];jo), 

Jl<...<j'n=l 

the expectation value of the number operator is easily 
calculated to be (n k ) n = |C n | 2 (n!) 2 |afc| 2 <S[ n ~ 1 , where 



s ( k n} = 



2M 

E 



(Cll) 



31<--<3'«=1 



If < \v k \ 2 < 1 — e, there exists a lower bound on the co- 
efficients K| 2 = M7\A ~ M 2 > 6 V fc. Then 



and 5*^™ 2 ' ) are related via 



M 2M-M (n _ 2) _ 
n - 1 
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In the proof of Thm. IC.2I we show that 

|A„| 2 |C„_i| 2 ((n-l)!) 2 



|Ar. 



|C„| 2 (n!) 2 



(C12) 



resulting in 



(nfc)n - (ifc)n-i > 

2M-(n-l) |A„| 2 



n - 1 



|C„| 2 (n!) 2 |a fe | 2 S 



2o(«-2) 



(C13) 



For n = N + A and A £ [ — CatjCat] the normal distri- 
bution of the |A„| 2 implies |A„| 2 /|A„_i| 2 = exp[(2A — 
l)/(2o&)] < e. Hence, ^ ^--(n-i) _ JML > ^ 6 > 

e 2Af-""-l) > 3 2M-(n-l) - FOT M = ~ X )> thiS iS 

equivalent to |a/c| 2 > 3/(2g — 1), which can be achieved 
for q 3> 1. The last condition is satisfied, as we are con- 
sidering dilute systems, where M ^> N. Thus, (rifc)„ — 
(nfc)n-i > 0- implying £ fc K| 2 ((n fc )„ - (n fe ) n _i) < 1, as 
\vk\ 2 < 1- Using (|C8|) and a telescope sum, we conclude 
that 



Ti < 



< 



E a* 



iv + A| 2 |A| 



£(|A 



2 |A| 



16 



v/2T 



(C14) 



Next, we derive the bound for the operator W. Its 
expectation value is given by 



(W) n = \C n \ 2 (n!f2j2 



<Xk\ 



2M 

E 



31<-<3 n -l 



(C15) 

For b6 [N — A, JV + A], we use the same argumentation 
we have used for bounding {rik)n — (^fe)n— l; 

to obtain 



(W) n - (W) n -i < 2. 



(C16) 



Further, (n k ) n = (pIpI_ m + h.c.) n /2 + {n k n k+M ) n due 
to the symmetry a k — a k +M- Hence, (W} n < 2n. Thus, 
up to a factor of 2 we obtain the same bound as for in- 
putting all the pieces together we find that 



\(H({v k })) var - (H({v k })) N \ < 

72 / 
2(1 - e) + —=a N + 12N ( 1 - Eif(a N /V2) 



(C17) 

In the limit of large x, the error function Erf (x/y/2~n) can 
be approximated by the following formula: 

1 - Erf(x/V2) = 2 eXp[ 7^L 2/2] (x- 1 - x- 3 + . . .). (C18) 
v 2tt 



As on = 0(yN), we conclude 

12JV(1 - Etf(af}/V2)) < 24a R exp[-cr 2 ? /2]/(eV2^ : ) -> 

for N 1. A straightforward calculation results in 
(H({v k })) var = — 47Ve, leading immediately to the state- 
ment of the theorem. □ 

APPENDIX D: LEMMA FOR THE PROOF OF 
THM. El 

Lemma D.l. Every pure separable state p € S{An) 
fulfills 53fcf_i \ {P k Pi) p\ ^ Af/2, an d this bound is tight. 
Proof. Using Lemma lA.2| we obtain 

M M 

E K p * tp <>i = E \ p M p -K-i - Pk-iP-ui (di) 



k,l=l 



k,l=l 



where P = P 2 = P^ and tr[P] = JV. Using the triangle- 
inequality we get 



a/ 



E k^>i 

k,l=l 

< \ E p «i 2 + i p -^i 2 + i p ^i 2 + i p -^i 2 ) 

= itr[P 2 ] = iV/2. 

In the last step we have used the property that the sum 
of the squares of a normal matrix is equal to the sum 
of squares of its eigenvalues. Taking the square root we 
obtain the bound of our claim. 

The bound is tight, as P — l 2jv implies J^ki 1(^^)1 = 



N/2 which is obtained for |$) = l\f=i 



APPENDIX E: PROOF OF LEMMA IHL51 



,t„t 



□ 



a'AO) and consider the subspace 



* 3 



Proof. Let \i,j) 

spanned by the states {\k, — k), \l, — I), \k, I), \k, — I), | 
k, l), | — k, — I}}. In this basis, the two-particle RDO p^ 
of I^Jgcs) 1S °^ ^ ne f° rm 



2ai 



a\ a-i 
a 2 a\ 
1 4 



(El) 



where a 1 = (M — 1)/(N — l),a 2 = (M — N)/(N — 1). 
The witness operator of Thm. IV. II has a negative 

expectation value on p^ N \ hence the state is paired in 
these modes. 

For solving the entanglement question we will use the fol- 
lowing theorem [25j applicable to mixed fermionic states 
of two particles each living on a single-particle Hilbert 
space of dimension four: 
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Theorem E.l. Let the mixed state acting on A4 have a 
spectral decomposition p — ^I=i l^i)(^i|> where r is the 
rank of p, and the eigenvectors belonging to nonzero 
eigenvalues Xi are normalized as ($>i\tyj) — XiStj. Let 
= ~^2 a b Waba^allQ) in some basis, and define the 
complex symmetric r x r matrix C by 



where 



-1^4 



5+a. 



■ I ip- ) and 



±k,±i) = 



' ~~ " 5+a + 

\l, —I)) and a+ 
7 2 = l/(5 + a + : 



±k,±l). Here |V±> = ^{\k,-k) ± 
(2M - N - 1)/(N - 1). Defining 



one obtains 



abed 



(E2) 



which can be represented using a unitary matrix as C — 
UCdU T , with Cd — diag[ci , . . . , cy] diagonal and \ci\ > 
I C2 1 > • ■ ■ > I Or- j - The state has Slater number 1 if and 
only if 



C = l 2 



( a+ \ 
0-10 
-1 
-1 
1 

V 010/ 



(E4) 



i=2 



The spectral decomposition of p% is given by 



■|*-><*-| + 
*m}(*m| + \^ k -i)(^k-i\ + 

*-«)(*-w| + \$-k-l)(*-k-l\ 



(E3) 



with spectrum spec(C) = 7 {a+, 1, 1, — 1, — 1, — 1}. For 



M < N the state I* 



(JV) 
BCSI 



is separable, so we can take 



M > N. Hence, a+7 2 is the eigenvalue with biggest 
absolute value. According to Thm. IE. 11 the reduced 
state in the subspace of the four modes is entangled iff 
|ci| < 7= 2 I c * I - F° r our example, this holds iff M > 
3N-2. □ 
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